Gravity Effects on Neutrino Masses in Split 

Supersymmetry 



Marco Aurelio Diaz, Benjamin Koch, and Boris Panes 
Departamento de Fisica, Universidad Catdlica de Chile 
Avenida Vicuna Mackenna 4860, Santiago, Chile 
(Dated: June 26, 2009) 

The mass differences and mixing angles of neutrinos can neither be explained by R-Parity vio- 
lating split supersymmetry nor by flavor blind quantum gravity alone. It is shown that combining 
both effects leads, within the allowed parameter range, to good agreement with the experimental re- 
sults. The atmospheric mass is generated by supersymmetry through mixing between neutrinos and 
neutralinos, while the solar mass is generated by gravity through flavor blind dimension five opera- 
tors. Maximal atmospheric mixing forces the tangent squared of the solar angle to be equal to 1/2. 
The scale of the quantum gravity operator is predicted within a 5% error, implying that the reduced 
Planck scale should lie around the GUT scale. In this way, the model is very predictive and can be 
tested at future experiments. 



I. INTRODUCTION 



The existence of neutrino masses and mixing angles is experimentally well confirmed and therefore 
the theoretical understanding and description of those quantities is one of the most urgent issues for particle 
physics I ill- There are two frequently studied theoretical extensions of the standard model of particle physics 
which also were expected to explain the origin and the shape of the small but non zero neutrino mass matrix. 

On the one hand there is Supersymmetry. Although there is no experimental evidence for supersym- 
metry, it is often invoked that the unification of gauge coupling constants at the GUT scale is an indirect 
indication for supersymmetry. It has been pointed out that this gauge unification is achieved even if all 
scalar superpartners of quarks and leptons are very heavy, introducing the scenario called Split Supersym- 
metry (SS) 113, yO. The original SS scenario includes R-Parity conservation, which guarantees the stability 
of the LSP and thus a Dark Matter candidate, and only one light Higgs doublet, which behaves like the SM 
Higgs. Despite the fact that models with R-Parity Violation (RpV) loose the LSP as a dark matter candidate, 
they are studied because they provide a compelling mechanism for the generation of neutrino masses 141]. 
Nevertheless, SS models with RpV are incapable to produce the necessary neutrino masses, even in next to 
leading order calculations ISilaLZl]. 
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On the other hand there is the possible existence of non-renormalizable gravitational interactions. Those 
interactions could have an influence on the neutrino sector of the standard model |^ ^ [l^ [l^ 13]. 
Although the standard Planck scale Mp = 1.2 x 10^^ GeV generates a solar mass which is too small to fit 
the experimental evidence, a lowered Planck scale Mj might in principle do the job. However, gravitational 
interactions are expected to be "flavor blind" and it has been shown that therefore a pure ly g ravitational 
neutrino sector is also incapable to explain the existence of three different neutrino masses 11411 . 

In this paper we study the combined effect of the supersymmetrical and gravitational neutrino sector. 
We find that the combination of both effects can explain all present neutrino data. But this kind of model 
even leads to predictions for the solar neutrino mixing angle sin^ 9 sol and for the scale of the gravitational 
contribution. 



II. THE NEUTRINO MASS MATRIX IN R-PARITY VIOLATING SPLIT SUSY 



In Split SUSY all scalars are very heavy, except for one Higgs doublet [2]. Since we are interested 
in the possibility to describe the neutrino masses in Split-SUSY, we need to consider R-Parity violating 
interactions jl5]. Knowing all relevant trilinear R-parity violating interactions one can calculate the neu- 
tralino/neutrino mass matrix. Integrating out the neutralinos from this matrix, one finds that the neutrino 



mass matrix in flavor space is given by 



M 



eff 



4 det M^o 



Nh~gl + M2~g'i) 



Af A1A2 A1A3 
A2A1 A2 A2A3 
A3A1 A3A2 A3 



where the determinant of the neutralino mass matrix is: 



(1) 



(2) 



Here, v is the vacuum expectation value of the light Higgs field, M\, M2 are the gaugino masses. Further, 
gu,d and g'^ ^ are the trilinear couplings between the Higgs boson, the gauginos, and the higgsinos. The 
parameters A^ = Oj/i + ej la] are related to the traditional BRpV parameters Aj [16] by Aj = XiVd- The 
are the parameters that mix higgsinos with leptons, and Oj are dimensionless parameters that mix gauginos 
with leptons. Finally, jj, is the higgsino mass. 

This effective neutrino mass matrix liAf,^-^ has only one eigenvalue different from zero. As in the case 
of R-parity violation in the MSSM with bilinear terms, only one neutrino is massive. As it is explained in 
the literature isi ^, in SS it is not possible to explain the neutrino masses and mixing using bilinear terms 
only. Further allowing for trilinear couplings makes it in principle possible to obtain solar neutrino masses 
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and mixing. However this possibihty was not discussed here, because it only works for a very special 
choice of parameters with an undesired hierarchical structure among the trilinear couplings 1^. Without 
this hierarchical structure, the trilinear couplings become irrelevant because they contribute through loops 
of sparticles that have a mass of the order of the split supersymmetric scale fh, and therefore, decoupled. 



III. THE NEUTRINO MASS MATRIX FROM LOW SCALE GRAVITY EFFECTS 



It is widely assumed that the unknown quantum gravity Lagrangian can be expanded to low energies. In 
flavor space this expansion can give rise to a non-renormalizable term in the Lagrangian of the type 



(3) 



Here, ipi and (p are the lepton and Higgs fields respectively. In an idealized model the flavor mixing coef- 



ficients Xij can all be taken to be of order one 
expectation value, a neutrino mass is generated 



m 



1411 . When the Higgs acquires a vacuum 



m,, ~ 



0{1] 



Mp 



0(W~n eV 



(4) 



with V the electroweak vacuum expectation value. This type of contribution to the neutrino mass matrix 



has also been explored in 111711 . Since this neutrino mass term is too small, we look into the possibility that 
the true Planck scale Mf is actually much lower than Mp. Such a lowered Planck scale M f is an intrinsic 
prediction of models with compact extra dimensions. In some of those models the experimentally allowed 
Planck scale Mf can be as low as one TeV HQHl. However, a TeV gravity scale in operators hke in 
eq. (O meets strong constraints from precision measurements such as ^ — > 67 ^J]. This can be met by 



imposing additional symmetries or by admitting that Mf might not be so extremely small 



1 TeV <^Mf < Mp 



(5) 



Several papers 



Q[n|,[i2, 



1411 discuss contributions of an exact "bhndness" model, where the part of 



the neutrino mass matrix coming from gravitational effects can be parameterized in flavor space as 



AM^ 



1 1 1 
1 1 1 
1 1 1 



(6) 



where 



/ig ~ ~ 0(10-^) eV 



(7) 
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It should be to noticed that such an exact "bUndness" model does not imply base independence in flavor 
space. This means that only when written in the flavor base (defined by a diagonal mass matrix of the 
charged Leptons), the matrix (O takes its symmetric form. Further, such a model can only give direct pre- 
dictions for the UPMNS matrix ( [TT] ) if the relation of flavor basis to the mass basis is defined for the charged 
leptons as well. Like in the standard model this relation is implicitly assumed for the above contribution by 
taking both basis to coincide 11211 . 

In order to explain the solar mass scale the fundamental scale has to be M f ~ 0(10^^) GeV, which is 
in good agreement with the inequality in eq. (O. However, the matrix in eq. Q has only one eigenvalue 
different from zero. Therefore such a scenario can not account for both solar and atmospheric neutrino mass 
splittings. 

IV. THE NEUTRINO MASS MATRIX FROM R-PARITY VIOLATING SPLIT SUSY AND LOW SCALE 

QUANTUM GRAVITY 

Now we will discuss the possibility that both, R-parity violating Split Susy and low scale quantum 
gravity operators like eq. ^ contribute to the Lagrangian. In this case, two terms can contribute to the 
neutrino mass matrix. Since the mass terms in eqs. ([T]) and (O are both formulated in flavor space one can 
write the total effective neutrino mass matrix as, 

up^ =Ay\^ +Hg , (8) 

where we defined the factor in front of the mass matrix in equation ^ as A. As defined above, the neutrino 
mass matrix does not contain CP violating phases. In this case, if j4 > and fXg > the three neutrino 
masses are 

ruui = 



m,, = i (a|A|2 + 3/.g) - y (^|A|2 + - 4Ai2g\v x A|2 (9) 



^ (a|A|2 + 3^3) +y (^|A|2 + 3/x,)' - 4Af,g\v X Xp 



where we have defined the auxiliary vector v = (1,1,1). In the approximation where the ^g term is 
subdominant, the squared mass differences are given by. 



Amf„, = {mi - mi) = fi'g ' + (10) 

Aml^ = {ml - ml) = AA^ + 2Af,g{vX)^ + 0{fil) (11) 
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This shows that for a very small fig, the atmospheric mass scale is controlled by the parameter A and the 
solar mass scale is controlled by /x^. In this approximation, the eigenvectors are equal to, 



Vl 



V2 



V3 



vxX 

\v X A| 

Xx (v xX) , 
\Xx(vx A)| 



(12) 



and the matrix Upmns is formed with the eigenvectors in its columns. The neutrino mixing matrix is 
defined as 



UpMNS = U32U31U21 , 



(13) 



where U32 is a rotation matrix around the axis one (with cyclic permutations for the other two). 

At this point it is instructive to have a closer look at the matrix Upmns- This matrix appears in the 
Lagrangian of the leptonic charged current interactions such that. 



(14) 



where fermions are in the mass eigenstate basis. In the most general situation, both the charged lepton and 
neutrino mass matrices are non-diagonal in the interaction basis. 



jr3-{e'+,fi'+,f'+)Me 



I- 



{i?[',i?'2',u'3')M, 



(15) 



where the prime on the fermion fields denote the interaction basis. The two mass matrices are diagonaUzed 
by Ve and K, 



V^MiVi = diag(me, m^, m^) , V^M^V^ = diag(m,,i , , rriv^) 



(16) 



and the U pm ns matrix becomes. 



UpMNS = V}Vi,, 



(17) 



where we have assumed CP conservation. Due to the Majorana nature of neutrinos, Upmns depends on 
three mixing angles and three Majorana phases. We assume the latter to be zero. In general, a simultaneous 
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TABLE I: Experimental measurements for neutrino parameters 



Observable Oi 


Mean value Oi 


3ai variance 


Units 




2.35 X 10-3 


0.95 X 10-3 


eV2 




8.15 X IQ-'^ 


0.95 X 10-5 


eV2 


sin 6atm 


0.51 


0.17 




sin^ 6 sol 


0.305 


0.075 





TABLE II: Experimental upper bounds for neutrino parameters 



Observable Oi 


3ai upper bound 


Units 


sin Oj-eac 


0.047 
0.84 


eV 



diagonalization of M,y and Mi is necessary to find Upmns- In principle there might be also off-diagonal 
contributions to the matrix coming from the "flavor blind" gravitational terms in eq. ^ (R-parity vi- 
olating supersymmetric intercations do not modify Vi). However, those contributions (if they are present) 
can be depreciated in our model since they only can produce corrections to the diagonal mass matrix of 
the charged leptons of the order ^ig/rrie ~ IQ-^. This makes Vi equal to unity up to terms of the order 
of i^Lg/me- In order to have a predictive model we do not consider any other beyond the standard model 
contributions, like non-diagonal terms in the Yukawa couplings of the charged leptons. Therefore, one can 
identify Upmns with Vy. It is in this basis that we include the "flavor blind" contribution from gravity 
given in eq. Q. 



V. RESULTS AND PREDICTIONS 



In this section we present some numerical results, where we find eigenvalues and eigenvectors of the 
3x3 effective neutrino mass matrix using numerical methods. With them, we find neutrino mass differ- 
ences and mixing angles, and compare them with values from experimental measurements. The solar and 
atmospheric mass differences, as well as the solar and atmospheric mixing angles, have been measured in 
several experiments. We use the results of the combined analysis in ref. [22], given in Table I. In addition. 



upper bounds have been obtained for the reactor angle i22h . and for the neutrino less double beta decay mass 



parameter m^^. The value for the neutrino less double beta decay parameter is quite precisely determined 
in our setting. Our model allows values of 0.001 eV< < 0.005 eV which is two orders of magnitude 
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TABLE III: Predictions for neutrino parameters 



Observable 


Value 


Units 




2.39 X 10-3 


eV2 




7.74 X 10-5 


eV2 


sin datm 


0.596 




sin^ 0soi 


0.321 




sin ^reac 


0.023 




"1/3/3 


0.0039 


eV 



below the cuirent experimental bound. 

We scan the parameters space, varying randomly the parameters Aj, A, and fig, and calculate the good- 
ness of the model, represented by eq. ([8]l, with 

i 

where we have assigned a null mean value to the two parameters in Table II for which only upper bounds 
are known. Positive and negative solutions for A were found. A typical solution for negative A is 

Ai = 0.0148 GeV^ , A2 = 0.0822 GeV^ , A3 = -0.0712 GeV^ , 

^ = -4.10 eV/GeV^ , = 0.003 eV (19) 

with = 0.02. This shows that it is possible to find solutions for the combined model which are in (3a) 
agreement with every single neutrino observable. The prediction for the neutrino parameters in this case are 
summarized in Table III. 

Based on the same Scr-limits scan, we study how the parameters of this model are constrained by the 
experimental data. A key ingredient of this model is that it offers a possible way to measure ij,g for "flavor 
blind" gravity effects. Indeed, an interesting prediction is that the allowed values of the mass parameter /ig 
are strongly constrained. In Fig. [T] we show the frequency of occurrence of each value of fig among the 
selected models. This model predicts values centered around fig = 3 x 10^^ eV with a 5% error, as we can 
see from the green region defined by < 1- For comparison we show also the red region for < 2, and 
the blue region for < 3. According to eq. (|7]l, this implies a prediction for the true Planck scale given by 
Mf PS 2 X 10^^ GeV, which is remarkably similar to the GUT scale. 

The same scan is shown in Fig. |2] where we see the neutrino mixing angles in correlation with the 
gravitational neutrino mass parameter fig. Concentrating on the green triangles (x^ < 1), we see that the 
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FIG. 1 : Prediction for the mass parameter /ig among models that satisfy all the experimental constraints. 
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FIG. 2: Sections of the parameter space showing sin^ 9soi and sin^ Oatm as a function of /ig, with < Ij 2, 3, for 
green triangles, red squares, and blue circles respectively. 



values of sin^ 9 sol (left frame) are centered around 0.33, which is in the larger side of the experimentally 
allowed window. On the other hand, the values of sin^ Oatm are centered around one half, which is right in 
the center of the experimental window. For comparison we show also points of parameter space less favored 
by experimental data, with red squares corresponding to < 2, and blue circles to < 3. 
Since the reactor angle satisfy. 



sm 



(^^3,1)^ 



< 0.047 



(20) 



where v^^i is the first component of the third eigenvector in eq. ((12)) . and the quoted upper bound corresponds 
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to the one given in Table II, we need ^ A2 + A|. Neglecting the value of Ai in front of A2 and A3 we 
obtain, 



2 



V3,lJ (A3-A2)2 2 



The numerical values for each parameter shown in eq. (I2TI ) are obtained in the following scenario. Exper- 
imental results indicate sin^ 9atm = 0.51 ± 0.17, where again we indicate the 3a error. Maximal mixing 
is satisfied if we take A| = Ag. Choosing the sign A3 = — A2 leads to the prediction for the solar an- 
gle indicated in the previous equation, sin^ Ogoi = 1/3, which nicely agrees with the experimental result 
sin2 Osoi = 0.305 ± 0.075. 




FIG. 3: Frequency of occurrence for the parameters Ai among the models compatible with experimental results. 



As we will see now, this prediction is confirmed by the scan and the numerical diagonalization of the 
neutrino mass matrix. In Fig. [3] we plot the frequency of occurrence for each of the parameters Aj among 
the models in our scan which satisfy the experimental constraints detailed in Tables I and II. We see in 
the left frame that Ai is typically an order of magnitude smaller than the other two parameters A2 and A3, 
consistent with requirements from the reactor angle. Due to this hierarchy in the parameters, our model 



has some interesting common features with the tri-bi-maximal mixing model discussed in [121]. In the 
right frame we see that models consistent with experiments need A3 = — A2, leading to the prediction 
sin^ 9 sol = 1/3 as explained earlier. The prediction tan^ ^23 = 1 ^ tan^ Ogoi = 1/2, is observed also in 
our scan. In Fig. |4] the relation between the neutrino mixing parameters sin^ 9atm and sin^ 9 sol is shown. 
From this plot one sees that the model can only deliver a good agreement with all experimental bounds, 
if 0.313 < sin^ 9soi < 0.342. This very strong constraint on sin^ 9soi (remember that the experimentally 
allowed region is 0.23 < sin^ 9soi < 0.38) gets even smaller if the atmospheric mixing parameter sin 9atm, 
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FIG. 4: Section of the parameter space showing sin^ datm versus sin^ 9soi, with < li 2, 3 for green triangles, red 
squares, and blue circles respectively. 



is taken at its central value of 0.51, as seen in Fig. 01 In this case one finds that the model predicts 0.325 < 
sin^ 9soi < 0.334. With those small theoretical uncertainties an improved measurement of sin^ 9soi would 
already allow to confirm the model prediction or otherwise rule out this model. 
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FIG. 5: Sections of the parameter space showing Ato^qj and Am^^^ as a function of fig with x'^ < 1- 



The neutrino mass differences Am^^^ and Am^^„ are shown in Fig. |5] The left frame shows that a 
growing ij,g gives a growing Am^^j, which can be understood by comparing to the approximation in eq. ([TOll . 
that is, the A term dominates over the /x^ term in eq. ([Hi. In the right frame we have the atmospheric mass 
difference Am^j„ as a function of the gravitational neutrino mass parameter fig, where there is no obvious 
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dependence because the atmospheric mass difference is dominated by the A term. In both mass differences 
though, we see that the predictions lie nicely at the center of the experimental window. 

It is also instructive to study how the scale of the supersymmetric term A and the scale of the gravitational 
term /ig are pinned down by the individual experimental constraints in Table I. Starting from the benchmark 
point in eq. ([T9l ) and varying with respect to A and fig we plot the allowed regions for every observable 
and the intersection of those regions. At the center of Fig. [6] we have the zone in the A- fig plane which is 




FIG. 6: Allowed regions in the positive A-fig plane. 



consistent with experiments, and formed by the intersection of several regions. Region I is a vertical stripe 
where mainly the A parameter is constrained by the Am^^^ data. Region II is a horizontal stripe where fig 
is constrained by the Am^^^ data. Note that this horizontal stripe turns into a vertical one after an eigenvalue 
crossing, although this last branch is not allowed by the solar angle data. Finally, region III indicates the 
constraint from sin^ 9 sol, whose allowed region is at the right of the diagonal line, which breaks at the top 
also due to an eigenvalue crossing. Constraints from the other neutrino observables are not relevant for this 
figure. 

The main features from Fig.[6]can be understood by looking at the approximations in eqs. (ITOl ) and (fTTI) . 
which show that for small values of fj,g, the atmospheric mass difference Am^^^ is dominated by A, while 
the solar mass difference Am^^^ is dominated by fig. One sees that in the vicinity of this benchmark point, 
the other observables play a minor role in pinning down the scale parameters A and fig. 
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VI. SUMMARY 



It is known that in Split Supersymmetry with R-Parity violation an atmospheric mass difference is gen- 
erated at tree level, but one-loop contributions are not enough to lift the symmetry of the effective neutrino 
mass matrix, thus not being able to generate a solar mass difference. This problem is not solved by adding 
the, in principle always present, "flavor blind" couplings from dimension five operators. The reason is 
that the Planck scale is too large to generate a solar mass difference large enough to be compatible with 
experiments. In this article we have shown that Split Supersymmetry with R-Parity violation, plus "fla- 
vor blind" gravity effects with a reduced Planck scale, present in models with compact extra dimensions, 
can be compatible with all data form neutrino experiments. The atmospheric mass difference is generated 
by supersymmetry with a mixing between neutrinos and neutralinos, while the solar mass difference is 
generated by the "flavor blind" gravitational effects. This model predicts a value for the gravitational term 
fig = 3x 10^^ lb 5% eV, which corresponds to a reduced Planck scale Mf 2 x 10^^ GeV. The fact that this 
reduced Planck scale is equal to GUT scale is a tantalizing result that may be related to gauge coupling unifi- 
cation of all four forces. In addition, the solar mixing angle is predicted to satisfy 0.313 < sin^ 0soi < 0.342. 
We show also that a maximal atmospheric mixing sin^ Oatm = 1/2 implies sin^ 6 sol = 1/3, which agrees 
with the implications from our parameter scan 0.325 < sin^ 0soi < 0.334, when we adopt the central value 
sin^ Oatni = 0.51. In this way, the model not only reproduce the experimental results but it is also predictive 
and, therefore, can be falsified by future experiments. 
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